Electroweak nuclear response in quasi-elastic regime 
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The availability of the double-differential charged-current neutrino cross section, measured by the 
MiniBooNE collaboration using a carbon target, allows for a systematic comparison of nuclear effects 
in quasi-elastic electron and neutrino scattering. The results of our study, based on the impulse 
approximation scheme and a state-of-the-art model of the nuclear spectral functions, suggest that 
the electron cross section and the flux averaged neutrino cross sections, corresponding to the same 
target and comparable kinematical conditions, cannot be described within the same theoretical 
approach using the value of the nucleon axial mass obtained from deuterium measurements. We 
analyze the assumptions underlying the treatment of electron scattering data, and argue that the 
description of neutrino data will require a new paradigm, suitable for application to processes in 
which the lepton kinematics is not fully determined. 
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The data set of Charged Current Quasi Elastic 
(CCQE) events recently released by the MiniBooNE col- 
laboration [l[ provides an unprecedented opportunity to 
carry out a systematic study of the double differential 
cross section of the process, 

v^.+^'^c ^ ^l- + x , (1) 

integrated over the neutrino flux. Comparison between 
the results of theoretical calculations and data may pro- 
vide valuable new information on nuclear effects, whose 
quantitative understanding is critical to the analysis of 
neutrino oscillation experiments, as well as on the ele- 
mentary interaction vertex. 

The charged current elastic neutrino-nucleon process 
is described in terms of three form factors. The vector 
form factors Fi{Q'^) and ^2(Q^) (Q^ = —q^, q being the 
four-momentum transfer) have been precisely measured, 
up to large values of Q^, in electron-proton and electron- 
deuteron scattering experiments (for a recent review, see, 
e.g., Ref.Q)- The Q^-dependence of the axial form factor 
Fa{Q^), whose value at = can be extracted from 
neutron /3-decay measurements, is generally assumed to 
be of dipole form and parametrized in terms of the so 
called axial mass Ma- 

Fa{Q^)^9a {l + Q^/Miy^ . (2) 

The world average of the measured values of the axial 
mass, mostly obtained using deuterium targets, turns 
out to be Ma = 1.03 ± 0.02 GeV d-H, while analyses 
carried out by the K2K Q and MiniBooNE coUab- 
orations using oxygen and carbon targets, respectively, 
yield Ma - 1.2^ 1.35 GeV. 

It would be tempting to interpret the large value of 
Ma reported by MiniBoonNE and K2K as an effective 
axial mass, modified by nuclear effects not included in 
the Fermi gas model employed in data analysis. However, 



most existing models of nuclear effects (for recent reviews 
see Ref.Q) fail to support this explanation, suggested by 
the authors of Ref. [7| , a prominent exception being the 
model of Ref. @. 

Obviously, a fully quantitative description of the 
electron-scattering cross section, driven by the known 
vector form factors, is a prerequisite for the under- 
standing of the axial vector contribution to the CCQE 
neutrino-nucleus cross section. 

Over the past two decades, the availability of a large 
body of experimental data has triggered the development 
of advanced theoretical descriptions of the nuclear elec- 
tromagnetic response. The underlying scheme, based on 
nuclear many-body theory and realistic nuclear hamilto- 
nians, relies on the premises that i) the lepton kinematics 
is fully determined and ii) the elementary interaction ver- 
tex can be extracted from measured proton and deuteron 
cross sections. 

The above paradigm has been successfully applied to 
explain the electron-nucleus cross section in a variety of 
kinematical regimes (for a recent review of the quasi- 
elastic sector see Ref. However, in view of the 
uncertainties associated with the energy of the incoming 
beam, the identification of the reaction mechanisms and 
the determination of the interaction vertex, its extension 
to the case of neutrino scattering may not be straightfor- 
ward. 

In this work we compare theoretical results obtained 
from the approach described in Refs. [HI, [13 to the mea- 
sured CCQE cross sections of Ref.[l[, discuss the differ- 
ences involved in the analyses of electron and neutrino- 
nucleus scattering, and argue that modeling neutrino in- 
teractions may require a paradigm shift. 

Electron-nucleus scattering cross sections are usually 
analyzed at fixed beam energy, i^e, and electron scatter- 
ing angle 0e, as a function of the electron energy loss w. 
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As an example, Fig. [T] shows the double differential cross 
section of the process 



X 



(3) 



at = 730 MeV and 9^ = 37°, measured at MIT- 
Bates The peak corresponding to quasi-elastic (QE) 
scattering and the bump at larger w, associated with ex- 
citation of the A-resonance, are clearly visible and well 
separated. Note that the three-momentum transfer |q| 
turns out to be nearly constant, its variation over the 
range shown in the figure being < 5%. As a consequence, 
the cross section of Fig H] can be readily related to the lin- 
ear response of the target nucleus to a probe delivering 
momentum q and energy w, defined as 



\{n 



2(5(6^ + So - S„) . (4) 



In the above equation, |0) and |7i) represent the target 
initial and final states, having energy Eq and i?„, respec- 
tively, while Ok+q and ak are the nucleon creation and 
annihilation operators. 

In addition, the magnitude of the momentum transfer, 
|q| ~ 450 MeV, is large enough to make the impulse 
approximation (lA) scheme, in which Eq.(j4]) reduces to 

m 



5/A(q,w)= j d^k dEPhik,E)Ppik + ci,u;-E) 



(5) 



safely applicable [lal- In Eq.(l5|), Ph{ k , E ) and 
Pp(k -f q , a; — E) are the spectral functions de- 
scribing the energy and momentum distributions of the 
struck nucleon in the initial (hole) and final (particle) 
states, respectively. 
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FIG. 1: (color online). Inchisive electron-carbon cross section 
at beam energy Ee = 730 MeV and electron scattering angle 
= 37°, plotted as a function of the energy loss ui. The data 
points are taken from Ref. p^ . 

The solid line of Fig. [T] represents the results of a the- 
oretical calculation of the QE contribution, carried out 
within the lA using the hole spectral function of Ref. [l^ 
and the recent paramctrization of the vector form factors 



of Ref. Final state interactions (FSI) between the 
struck nucleon and the recoiling spectator system [Tl| . 
whose main effect is a ^ 10 MeV shift of the QE peak, 
have been also taken into account. 

It is apparent that height, position and width of the 
QE peak, mostly driven by the energy and momentum 
dependence of the hole spectral function, are well repro- 
duced. 
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FIG. 2: (color online). Flux integrated double differential 
CCQE cross section measured by the MiniBooNE collabora- 
tion [l[ , shown as a function of kinetic energy of the outgoing 
muon. The upper and lower panels correspond to to different 
values of the muon scattering angle. Theoretical results have 
been obtained using the same spectral functions and vector 
form factors employed in the calculation of the electron scat- 
tering cross section of Fig. [1] and a dipole parametrizaition 
of the axial form factor with Ma ~ 1-03 MeV. 

Applying the same scheme employed to obtain the solid 
line of Fig. [T] to neutrino scattering one gets the re- 
sults shown in Fig. [21 The data points represent the 
double differential CCQE cross section integrated over 
the MiniBooNE neutrino fiux, whose average energy is 
{ Eu ) = 788 MeV, plotted as a function of the kinetic 
energy of the outgoing muon at different values of the 
muon scattering angle. The solid lines show the results 
(integrated over the cos 0^ bins) obtained using the same 
spectral functions and vector form factors employed in 
the calculation of the electron scattering cross section of 
Fig. [Tl and a dipole parametrization of the axial form 
factor with Ma = 1-03 MeV. Comparison of Figs. [T]and 
[5] indicates that the electron and neutrino cross sections 
corresponding to the same target and comparable kine- 
matical conditions (the position of the QE peak in Fig. 
[T] corresponds to kinetic energy of the scattered electron 
^ 610 MeV) cannot be explained using the same theo- 
retical approach using the value of the axial mass result- 
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ing from deuterium measurements. In this instance, the 
paradigm of electron scattering appears to conspicuously 
fail. 

Note that the comparative analysis of electron and 
neutrino data, based on double differential cross sections 
depending on measured kinematical variables, is made 
possible for the first time by the availability of the data 
set of Ref. 

The results of a global comparison between the Mini- 
BooNE data and the calculated cross sections show that 
theory sizably underestimates the measured cross section 
for any values of muon energy and scattering angle. 

It has to be emphasized that the above conclusion, 
while being based on a calculation carried out within the 
scheme of Refs. [HI is largely model independent. 
Theoretical approaches providing a quantitative descrip- 
tion of the electron-nucleus cross section in the QE chan- 
nel, are bound to predict CCQE neutrino-nucleus cross 
sections significantly below the MiniBooNE data if the 
value of the axial mass is set to 1.03 GeV. 
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FIG. 3: (color online). Flux integrated muon kinetic en- 
ergy spectrum. The dot-dash, solid and dashed lines have 
been obtained setting the value of the axial mass to Ma = 
1.03, 1.35 and 1.6 GeV, respectively. The data are taken from 
Ref.tl . 

In spite of the fact that the large value of Ma reported 
by the MiniBooNE collaboration was obtained from a 
shape analysis of the Q^-distribution, the effect of the 
axial mass on the CCQE cross section can be best an- 
alyzed studying the flux integrated muon kinetic energy 
spectrum and angular distribution, obtained integrating 
the double differential cross section over cos0^ and T^, 
respectively. These quantities only depend on the mea- 
sured muon kinematical variables, thus being unaffected 
by the assumptions associated with the reconstruction of 
the incoming neutrino energy E^, [igI . [17| , entering the 
definition of the reconstructed Q^. 

Figure [3] shows a comparison between the MiniBooNE 
flux integrated muon kinetic energy spectrum and the re- 
sults of our calculations, corresponding to three different 
values of Ma- The behavior of the curve corresponding 
to Ma — 1.03 GeV is consistent with that shown in Fig. 
[21 as the data turns out to be largely underestimated. 



Increasing Ma to 1.35 GeV, the value resulting from 
the MiniBooNE analysis of Ref.[l|, while improving the 
agreement between theory and experiment, still does not 
lead to reproduce the data at < 1 GeV. The dot-dash 
curve has been obtained using the value AIa = 1.6 GeV, 
yielding the best x^-fit to the MiniBooNE fiux integrated 
Q^-distribution within the our approach. 
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FIG. 4: (color online). Flux integrated muon angular dis- 
tribution. The dot-dash, solid and dashed lines have been 
obtained setting the value of the axial mass to Ma = 
1.03, 1.35 and 1.6 GeV, respectively. The data are taken 
from Ref. 01 . 

The M^-dependence of the flux integrated muon angu- 
lar distribution, is shown in Fig. |4l together withe data 
from Ref. [l|. The overall picture is clearly the same as 
in Fig. [21 

In Fig. [SI we compare the results of our calculations 
to the MiniBooNE fiux unfolded total cross section. It 
is apparent that in this case using Ma = 1.6 GeV leads 
to overestimating the data in the region of high energy 
{E^ > 800 MeV), where the choice Ma = 1.35 GeV 
provides a better fit. The different pattern emerging from 
Fig. [51 compared to Figs. [3] and [H clearly points to 
the uncertainty associated with the interpretation of fiux 
averaged and fiux unfolded data. 

The results of our work indicate that the theoretical 
approach based on the lA and realistic spectral func- 
tions, successfully applied to QE electron scattering, fails 
to reproduce the CCQE neutrino-nucleus cross section, 
unless the value of the nucleoli axial mass resulting from 
deutcron measurements is modified. In addition, the pos- 
sibility of interpreting the large Ma resulting from the 
K2K and MiniBooNE analyses as an effective axial mass, 
modified by nuclear effects beyond the Fermi gas model, 
appears to be ruled out [l3|. This statement should be 
regarded as largely model independent, as calculations 
carried out using different descriptions of nuclear dynam- 
ics yield similar results [isj . 

A different scenario is suggested by the results of 
Ref. 0, whose authors obtain a quantitative account of 
the MiniBooNE fiux unfolded total cross section without 
increasing Ma- Within the model of Ref. Q, the main 
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FIG. 5: (color online). Flux unfolded total CCQE cross sec- 
tion, as a function of neutrino energy. The dot-dash, solid 
and dashed lines have been obtained setting the value of the 
axial mass to Ma ~ 1.03, 1.35 and 1.6 GeV, respectively. 
The data are taken from Ref. [l|] . 

mechanism responsible for the enhancement that brings 
the theoretical cross section into agreement vifith the data 
is multi nucleon knock out, leading to two particle- tvifo 
hole (2p-2h) nuclear final states. 

Within the approach employed in our work, the oc- 
currence of 2p-2h final states is described by the contin- 
uum part of the spectral function, arising from nucleon- 
nucleon correlations It gives rise to the tail extend- 
ing at large w, clearly visible in Fig. [TJ However, its con- 
tribution turns out to be quite small (less than 10% of the 
integrated spectrum of Fig. [T]) . According to the philos- 
ophy outlined in this paper, in order to firmly establish 
the role of multi-nuclcon knock out in CCQE neutrino 
interactions the model of Ref. @ should be thoroughly 
tested against electron scattering data. 

In our opinion, the available theoretical and experi- 
mental information suggests that the main difference in- 



volved in the analysis of neutrino-nucleus scattering, as 
compared to electron-nucleus scattering, lies in the fiux 
average. 

Unlike the electron cross section shown in Fig. [TJ the 
flux averaged CCQE neutrino cross section at fixed en- 
ergy and scattering angle of the outgoing lepton picks up 
contributions from different kinematical regions, where 
different reaction mechanisms dominate. As a conse- 
quence, it cannot be described according to the paradigm 
successfully applied to electron scattering, based on the 
tenet that the lepton kinematics is fully determined. 

A new paradigm, suitable for the studies of neutrino 
interactions, should be based on a more flexible model of 
nuclear effects, providing a consistent description of the 
broad kinematical range corresponding to the relevant 
neutrino energies. 

Besides single- and multi- nucleon knock out, such a 
model should include the contributions of processes in- 
volving the nuclear two-body currents, which are known 
to provide a significant enhancement of the electromag- 
netic nuclear response in the transverse channel [l^. In 
addition, the occurrence of processes leading to pion pro- 
duction and excitation of nucleon resonances should also 
be taken into account. 

A great deal of information could be obtained apply- 
ing the new paradigm to the analysis of inclusive data, 
preferably, although not necessarily, through direct im- 
plementation in the Monte Carlo simulation codes. This 
kind of analysis may help to reconcile the different val- 
ues of the axial mass obtained from different experiments, 
and would be unaffected by the problem of the possible 
misidentification of CCQE events, recently discussed in 
Ref. 

The authors are indebted to R. Schiavilla for a critical 
reading of this manuscript. Useful discussions with G.T. 
Garvey are also gratefully acknowledged. 
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